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BRANCHING  PROCESSES 
by 

T.  E.  Harris 


X.  Sirnmarv.  This  paper  is  concerned  with  a  simple  mathematical  model 
for  a  branching  stochastic  process.  Using  the  language  of  family  tree* 
we  may  illustrate  the  process  as  follows.  The  probability  that  a  man 
has  exactly  r  sons  is  pr,  r  *  0,1,2, ...,  Each  of  his  sons  (who  to¬ 
gether  make  up  the  first  generation)  has  the  same  probabilities  of  having 
a  given  number  of  sons  of  his  own;  the  second  generation  have  again  the 
sane  probabilities,  and  so  on.  Let  zn  be  the  number  of  individuals  in 

• 

the  nth  generation.  Ve  study  the  probability  distribution  ofA«^,  €foa* 
previbusv  results  are  given  in  section  2;  these  include  procedures  for 
computing  moments  of  zR,  and  a  criterion  for  when  the  family  has  pro¬ 
bability  1  of  dying  out.  In  sections  3  and  4  the  case  is  considered 
where  the  family  has  a  nonzero  chance  of  surviving  indefinitely.  In 
this  case  the  random  variables  zn/Ezn  converge  in  probability  to  a  ran¬ 
dom  variable  w  with  cumulative  distribution  G(u).  It  is  shown  that  G(u) 
is  absolutely  continuous  f o  •  u^O,  Results  of  a  Tauberian  character 
are  given  for  the  behavior  of  G(u)  as  u  — ►  0  and  u  — ►■ao.  In  section  5 
some  examples  are  given  where  G(u)  can  be  found  explicitly;  G(u)  is 
computed  numerically  for  the  case  p^  =  0,4,  =  0,6,  In  section  6 

families  with  probability  1  of  extinction  ire  considered.  A  method  is 
given  for  obtaining  in  certain  cases  an  expansion  for  the  moment-gener¬ 
ating  function  of  the  number  of  generations  before  extinction  occurs. 
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Xk  section  7  maximum  likelihood  estimates  are  obtained  for  the  and 

r 

f 

for  the  expectation  Ez^;  consistency  in  a  certain  3ense  is  proved.  In 
section  8  a  brief  discussion  i3  given  of  the  relation  between  two  types 

of  mathematical  models  for  branching  processes. 

% 

2.  Introduction.  By  a  branching  stochastic  process  is  meant  a  phenome¬ 
non  of  the  following  general  type:  each  of  an  initial  aggregate  of  ob¬ 
jects  can  give  rise  to  more  objects  of  the  same  or  different  types;  the 

objects  produced  can  then  produce  more,  and  the  system  develops,  subject 

■ 

« 

to  certain  probability  laws.  Examples  are  the  development  of  human  or 
animal  populations,  propagation  of  genes,  and  nuclear  chain  reactions. 
The  mathematical  model  dealt  with  in  this  paper  may  be  thought  of  as 
representing  the  generation-by-generation  growth  of  a  family,  the  funda- 
mental  random  variable  being  the  number  of  individuals  in  the  n1'  gener¬ 
ation.  Under  certain  conditions,  however,  this  model  may  describe  the 
size  of  a  family  at  a  sequence  of  points  in  time.  This  question  will 
be  touched  on  in  section  8, 

Definition  2.1.  The  random  variables  zR,  n  *  0,1,2,...,  will  be 

said  to  represent  a  simple  discrete  branching  process  provided:  z  =  l* 

co  0 

P(zl=r)  -  Pr,  r  r  0,1,2,...,  with  >  p  «  1;  the  conditional  distri- 

r*0 

bution  of  given  zn=r,  is  that  of  the  sum  of  r  independent  random 

variables,  each  having  the  same  distribution  as  z.. 

1  00  2 

Assumptions.  Throughout  this  prper  we  assume  that  ^  r  p  <  oo, 

r*o  r 

that  at  least  two  of  the  pr  are  positive,  and  that  P0+P^  <  1* 

Definitions  2.2,  Let  x  =  Ez^  -  rpr,  tf2  =  VarCz^)  * 

5  r  P  -x  •  Let  f (s)  =  ^  p rsr  be  the  generating  function  of  z. 

r=o  ^ 
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(s  denotes  a  complex  variable.)  Let  p  *  P(z  *r)  and  f  (s) 

nr  n  n 

oo 

/  p  s1 ;  of  course  p.  *  p  and  f  (s)  •  s.  The  assumptions  given 
r—o  nr  xr  r  o 

above  insure  that  the  first  and  second  derivatives  f*(s)  and  f  (s)  art 
continuous  in  the  set  consisting  of  the  interior  of  the  unit  circle  and 
the  point  s  =  1 ;  thus  derivative  notations  such  as  fM(l)  are  used  even 
though  f(s)  may  not  be  analytic  at  s  =  1.  It  will  be  seen  shortly  that 
a  similar  remark  applies  to  the  functions  fn(s)  and  certain  functions 
to  be  introduced  later. 

In  the  remainder  of  this  section  iae  shall  summarize  certain  results} 
most  of  them  are  contained  implicitly  or  explicitly  in  works  by  Fisher  [l], 
Lotka  [2]  ,  Ftef fensen  [3]  ,  Ulam  and  Hawkins  [4J  ,  Kolmogoroff  [5]  , 
Kolmogoroff  and  Dmitriev  [6],  and  Yaglom  [7];  some  of  these  references 
are  not  widely  available. 

From  our  definition,  P(zn+1*k  zn«j)  ls  th®  coefficient  of  sk  In 
[f(s)]J.  Henoe  p^^  ^  is  the  coefficient  of  sk  in  ]>  Pnj  |f(s)] 

whence 


(2.1) 


fn/)!fnN- 


Letting  n  *  1,2,...,  successively,  it  follows  that  the  generating  function 
of  zn  is,  the  nth  functional  Iterate  of  f(s).  Hence 

(?.?)  fn+1(s)  =  r[fn(s)]  . 

Ve  note  that  f^(l)  =  Ezn,  f^'(l)  +  f^(l)  -  ^(l)]2  »Var(zn).  Differ- 
entiation  of  (2.1)  at  s  =  1  gives  f^U)  r  *  ;  another  differentiation 
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gives  s  f”(l)  [f^(l)]  2  ♦  f*(l)f”(l)  while  twofold  differentiation 

of  (2.2)  gives  f”x(l)  =  f"(l)f^(l)  f  [V(l)]2  f”(l);  these  two  expres- 

»t  ti 

sions  for  fn^x(l)  can  be  equated  and  solved  for  fR  (1),  provided 

x  =  f'(l)^  1.  Thus  the  mean  and  variance  oL  zn  are  given  by 

E*n  =  (  Ez1)n  =  xn;  Var(zn)  =  (X2xn(xn-1) .  x^l$  Var(zn)  =  ,  x  =  1. 

x2-x 

Higher  moments,  if  they  exist,  may  be  found  by  a  similar  process. 

Definition  2.3*  Denote  by  a  the  smallest  non-negative  real  root 
of  the  equation  t  -  f(t).  We  see  that  x:^l  implies  a  -  1  while  x>l 
implies  O^a  <1,  the  equality  a  =  0  holding  if  and  only  if  p^  =  0. 

In  no  case  can  the  half-open  interval  0^t<l  contain  more  than  one 
root.  It  is  readily  seen  that 

(2.3)  lim  pno  =  lim  fn(0)  =  a. 
n— *oo 

We  thus  have  the  well  known  results  the  number  a  is.  the  probability  of 
eventual  extinction  of  the  family.  The  relation  between  a  and  x  shows 
that  the  probability  of  extinction  is.  1  if  and  only  if  x  £  1. 

It  is  also  clear  that  0  t  <  1  implies  lim  fR(t)  r  a;  this,  to¬ 
gether  with  (2.3),  shows  that 

(2.4)  lim  pnr  =  0,  r  -  1, 2,.... 
n — Voo 

Relation  (2.4)  means  roughly  that  the  family  either  dies  out  or  gets 
very  large.  In  section  4  it  will  be  shown  that  (2.4)  holds  uniformly 
in  r. 

Definition  2.4.  The  random  variables  w  are  defined  by  s  z  /x11. 

■  n  n  n 
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Clearly  Ew  2 


2  1  and  Ew2  »  1  +  ( 1-  l  \  if  x=X  1. 

n  fcv  **v  T 


Suppose  n>  m.  Then  E(znzm)  *  £Z  pmrE^rzn 

r 

m  2  2 

x  Ez  .  Thus  E(w  w  )  =  Ew  ,  whence 
m  n  m  nr 


z  =r) 
m  ' 


K —  2  n-m 

H  Pmr*"  * 


(2.5) 


2  2  2 

E(w  -w  )  *  Ew  -  Ew_.#  n  >  m. 

n  m  n  mr  ^ 


By  virtue  of  (2.5)  we  obtain 

Theorem  2.1.  I£  x  1,  the  random  variables  wn  converge  in  mean 

square,  hence  in  probability,  to  a  random  variable  w. 

2  rr^ 

For  in  this  case  Fw* — VI  4  2- —  as  n  - V  Oo  and  (2.5)  shows  that 

n  x-x 

2 

E(wn-wm)  - >  0  as  n  and  m - Voo.  Theorem  2.1  is  then  a  consequence 

of  [8],  p.  38,  I. 

It  is  well  known  that  convergence  in  mean  square  implies 


Ew' 


n 


Ew2  and  F(wn~i)2 


4-  E(w-l)  whence  Ew, 


Ew. 


Thus  we  have 


(2.6) 


Ew  =  1,  Ew2  =  1  t  2^ —  . 


X  -X 


In  order  to  study  the  behavior  of  z„  for  large  n  when  x  >  1,  we  consider 
the  distribution  of  w. 

Definition  2.5.  G  (u)  =  P(w  ^u);  0  (s)  =  e(c  n  )  =  /  eSUdG  (u), 

n  n  n  \  ©—  ® 

Definition  2.6.  (Applicable  when  x  >*1).  G(u)  *  P(w^u); 
roo 

0(s)  =  E(ews)  =  /  esudG(u).  We  shall  refer  to  G'(u)  as  the  asymptotic 

J  o- 

dlstrlbution  branching  fron  f ( s ) . 
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The  moment-generating  functions  (a.g.f.’s)  $n(£)  an<l  0(s)  are 
defined  at  least  for  Pe(s)^_;0.  Unless  specifically  stated  otherwise 
we  shall  consider  them  only  in  that  domain. 

From  (2.2)  and  the  fact  that  #n(s)  =  f n [jeS/^X  J  it  follows  that 

0n+^(sy)  =  f  pn  ( 2  )J  •  Theorem  2.1  implies  that  if  x>i  G^(u) - V  G(u) 

and  0n(s) - K0(s)  for  Fe(s)  ^0.  Thus  the  m.£.f,  0(s)  satisfies  the 

functional  equation 

(2.7)  0(sx)  ■  f  0(s)]  ,  Pe(s)^0. 

Equation  (2.7),  which  of  course  is  applicable  only  when  x>l,  was  ob¬ 
tained  in  a  different  form  by  Ulam  and  Hawkins.  It  belongs  to  a  type 
usually  known  as  Koenig’s  equation,  after  the  nineteenth  century  mathe¬ 
matician  who  studied  it  in  connection  vith  functional  iteration,  and  is 
related  to  an  equation  studied  by  Abel.  Y.e  shall  make  some  use  of  the 
work  of  Koenig’s  later.  See  Haderaard  [9]  and  Koenig  [lo]  . 

We  note  that  Ew  <.  OO  if  and  only  if  Ez^<0O.  It  was  already 

pointed  out  that  Ew  =  1.  As  pointed  out  in  [4],  as  many  further  moments 
of  w  as  exist  may  be  found  by  successive  differentiation  of  (2.7)  at 
s  -  0  • 

Finally  we  note  that  Gn(0)  *  pnQ.  Hence  Urn  Gn(0)  =  a*  Thus 

G(0)  -  P(w-0)^a.  Y/e  show  later  that  G(0)  -  a.  Clearly  G(u)  =  0 
for  u  <0. 

In  sections  3  and  4  we  always  assume  x  >  1 . 

3.  Asymptotic  properties  of  the  moment- genera  tint.  function.  We  first 
show  that  (2.7)  uniquely  determines  the  distribution  of  w. 
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Specifically, 

Theorem  3.1.  Let  G^(u)  and  G2(u)  ~  distributions  with  ecual  first 
moments  and  finite  second  moments  shore  characteristic  functions  0^(it) 

and  02(it)  satisfy  (t  is^  real)  0r(itx)  *  f  ^0^ (i t)J  ,  r  *  1,2,  Then 
Gx(u)  =  G2(u). 

From  jb]  ,  p.  27,  ^(it)  -  0j(lt)  =  t2/&(t),  where  /?(t)  is  bounded 
as  t - ►  0.  From  (2.7),  ^(ltx)  -  02(itx)  =  |  f  ^(it)]  -  fj02(lt)]|^ 

x|0^(it)  -  02(it)j,  since  f  (s)  <x  when  s  1.  Hence  for 

t  ^0,  ^  *  j^?(t)  |  .  Thus  /?(t)  cannot  be  bounded  near  t  =•  0 


unless  it  is  identically  zero;  hence  0^ (it)  =  02(it). 

It  is  clear  that  the  requirement  that  0(s)  have  the  form 
Its  *0(»2)  between  two  rays  from  the  origin  is  sufficient  for  the 
uniqueness  in  that  domain  of  solutions  of  (2.7),  On  the  other  hand, 
continuous  solu  ons  can  be  constructed  at  will  if  the  existence  of  a 


derivative  n 


Bef  ore 


0  is  not  required. 

ng  further,  it  is  convenient  to  define  three  functions 


k(s),  'f'(s),  ana  H(u)  which  are  closely  related  to  f(s),  0(s),  and  G(u) 
respectively.  We  repeat  that  we  are  considering  only  the  case  x>l. 


See  definition  2.3  for  a. 

Definition  3.1.  Let  k(s)  = 


1-a 


^iJL.  Clearly  k(s)  is  a 


probability  generating  function  with  k(0)  *  0,  k*(l)  s  f*(l)  *  x, 


k  (1)  <Oq«  We  write  k(s)  = 


q  s  .  We  also  define  the  iterates 
r 


hn(s)  by  k0(£)  =  s,  lcnU(s)  *  k[kn(sj] 
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Definltlon  3.2.  Let  H(u)  be  the  asymptotic  distribution  branching 
from  k(s).  (See  Definition  2.6).  Let  s )  be  the  corresponding 

moment-generating  function.  We  know  then  that  ^(s)  and  k(s)  satisfy 

0.1)  'f(sx)  =  k  [T<a>]. 

In  view  of  the  uniqueness  theorem  we  have,  by  direct  substitution  in 

(3.1) ,  that  ^(s)  must  be  given  by 

(3.2)  Y(s)  r  6f(l-a)J-  a, 

I- a 

and  that  H(u)  must  be  given  by 

(3.3)  H(u)  *  c(xi^)-a  ,  u>0;  H(u)  *  0,  u<0. 

1-a 


We  shall  see  later  that  H(0)  =  0;  i.e.,  that  G(0)  *  a.  Therefore 
H(u)  is  the  conditional  distribution  of  (l-a)w,  given  that  w^O. 
Another  way  of  stating  this  is  as  folio? s: 

Theorem  3.2.  The  random  variable  w  is;  distributed  as  the  product 


of  two  Independent  random  variables  w^  *  w  ,  where  wQ  takes  the  values 

0  2Si  IT.  with  probabilities  a  and  1-a  respectively  while  w*  has  the 
asymptotic  distribution  branching  from  k(s) . 

For  it  is  directly  verifiable  that  ^(s)  is  the  m.g.f.  of  wq*  w  . 

In  Theorems  3.3  ana  3.4  we  consider  the  behavior  of  ^(s)  for  large 
| s |  •  To  make  for  smoother  reading  we  defer  the  proofs  till  section  9, 
where  somewhat  more  general  formulations  are  given.  In  section  4  the 
properti  es  of  t'(s)  are  interpreted  in  terms  of  G(u). 

i 

(See  defini- 


Definltlon  3.3.  Let  “/  »  logx  f  q~)  s  p^" 


Ml 


tions  2.3  and  3.1).  If  q.  =0  (i.e.,  p  “p.sO)  we  take  - 

1  o  x 


Oo 
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Theorem  3.3*  Suppose  -y  <(J5  •  Then  Pe(s)  ^  0  and  s  ^4  0, 


(3-4) 


y(s)  =  ♦  M0(s). 

lslY 


H(s)  ir  continuous  for  s^Oj  M(?)  and  MQ(s)  satisfy  respectively 
(3.5)  M(sx)  =  V(s)|  U0(s)  :Q  f  1A  |sj - hoo- 

\  I*'  / 

Pemarks.  (See  section  9  for  proof.)  (a).  Under  the  conditions 
of  the  theorem  M(s)  is  real  and  positive  when  s  is  real  and  negative, 
(b)  If  Ez*  <OOand  the  conditions  of  the  theorem  hold,  the  rth  deri¬ 
vative  of  ^(s)  satisfies 


(c)  If  V  =  bo  ,  'f'(s)  and  as  many  derivatives  as  exist  approach  0 
exponentially  as  |s| - h  Oo* 

YJe  now  consider  the  behavior  of  ^ (s)  on  the  positive  real  axis, 
provided  it  is  defined  there. 

Lemma  3.1.  Let  f (s)  be  analytic  in  the  circle  |s|<TO(.,OC>  1* 
Then  0(s)  and  n  s)  r,re  analytic  in  some  neighborhood  of  s  •  0. 

Ve  use  a  theorem  of  Poincare""  [u]  which  insures  that  there  is 
exactly  one  function  0(s)  analytic  near  s  *  0  with  0(0)  =  (0)  •  1 

and  satisfying  0(sx)  *  f[0(s)]  .  (Although  Poincare'* s  proof  is  for  the 
case  f(s)  rational,  it  applies  equally  well  here.)  The  circle  of  con¬ 
vergence  of  the  MacLaurin  series  for  0(e)  has  radius  where 

)  =  CL,  An  arpument  whose  details  are  given  in  [l2]  ,  p.21,  then 
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shows  that  0(f)  -  ?(s)  for  s  <1 1^,  and  Lemma  3.1  follows.  (The 
argument  is  necessary  to  rule  out  the  possibility  that  the  0n(s)  con¬ 
verge  to  ?(s)  for  Pe(s)£^0  hut  to  some  other  function  for  Fe(s)>0.) 
Clearly  0(s)  and  ^(s)  are  entire  if  and  only  if  f(s)  is  entire. 

Lemma  3.L  is  usefuL  for  actual  computation  of  G(u).  The  (non- 

2 

negative)  coefficients  cr  in  the  series  0(s)  =  1  +  s  4  c, s  ♦ . . .  can 

be  determined  by  differentiating  (>.7)  at  s  =  0.  The  series  can  be 
used  to  compute  values  of  the  characteristic  function  0(it)  on  some 
interval  t^^t^t^y,  where  tQ  is  a  small  real  number;  the  values  of 

0(it)  for  the  remaining  values  of  t  are  determined  by  (^. 7 ).  (Note 
that  the  real  and  imaginary  parts  of  0(it)  are  respectively  even  and 
odd.)  Then  the  usual  inversion  formula  is  used  to  obtain  G(u).  A 
numerical  example  of  this  procedure  is  worked  out  in  section  5. 

Definition  3.4.  The  number  ^  is  defined  by  /°=  log  d  if  f(s) 
is  a  polynomial  of  degree  d,  -  0O  otherwise. 

Theorem  3.4»  Let  f(s)  (and  hence  k(s))  be  a  polynomial  of  degree 
d .  Then  for  s  0 


lOF 


si 


L(s)  4  L0(s); 


L(s)  continuous  and  positive:  L(s)  and  Lq(s)  satisfy  respectively 

L(sx)  -  L(s);  .j^(s)  =  '  5  ^  ^  * 

The  proof  is  in  section  (Theorem  3.4  may  be  compared  with  a  more 

widely  applicable  but  less  precise  result  due  to  Shiah  [1.9]  ) . 
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Corollary.  If  f(s)  is.  a  polynomial  of  decree  d,  s )  is.  an 

•entire  function  of  order  and  type  C  where  C  *  Max*  L(s),  1  Ss.  s  f^x. 

An  explicit  determination  for  C  has  not  been  found.  An  approxi¬ 
mate  numerical  determination  is  not  difficult;  the  function 

L(s)  -  lim  tn4  can  be  determined  numerically  for  a  number 

n — >0n 

of  values  on  some  convenient  interval  s 

< 

value  approximated.  The  importance  of  C  will  be  indicated  in  the  con- 
Jecture  folloving  Theorem  4 .3.  V.'e  may  also  mention  that  the  quantity 
JHax  L(s)  -  Min  L(s)]  ,  lf-^sf-^x,  is  of  some  interest.  Some  numerical 
work  indicates  that  in  certain  cases  L(s)  is  at  least  approximately 
constant . 


s  z^sqx,  and  the  maximum 


4.  Some  properties  of  C(u).  Since  it  will  be  convenient  to  work  with 
H(u)  rather  than  G(u),  ve  state  the  content  of  Theorems  4.1,  4.2,  and 


4.3  in  terms  of  G(u):  G(u)  *  a  ♦ 


•u 


g(v)dv  for  u>0.  Tne  density 


/  k  C  r )  * 

p(u)  is  continuous  for  u  =^0.  If  Ez  <C  00  then  gv  '  (u)  is  continuous 

/  1 

for  u  0  provided  r  <C  V  *  k  -  1  and  is  continuous  for  u  =  0  provided 
r  <V  -1.  Near  u  -  0,  G(u),  provided  V  <C0O,  approximates,  in  a 
certain  mean  sense  made  clear  by  Theorem  4.2,  the  function 

a  f  Lkr.sJ -  Mfu(l-a)]  ,  where  for  convenience  ve  have  defined  M(u) 

r  dtV)  L  J 

for  positive  u  by  M(u)  =  M(-u).  It  is  then  shown  that  in  a  certain  sense 

C- £  Q* £ 

p(u)  roes,  to  ?ero  faster  than  exp(-u*  )  and  slover  than  exp(-u  ) 
where  £  is  any  positive  number,  Q  being  defined  in  Theorem  4*3.  A  con¬ 
jecture  is  riven  of  a  more  precise  result,  applicable  when  f(s)  is  a 
polynomial:  in  the  same  sense  g(u)  goes  to  zero  (more, less)  rapidly  than 

(exp[-(A*-  £  )u<5]  ,  exp  ^  ( A  -fEju^J),  where  A  is  defined  in  the  conjecture 
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Deflnltlon  4.1.  Let  H  (u)  =  h(u). 

Theorem  4*1*  H(u)  ijs  absolutely  continuous.  Theorem  3.3  shows 
that  H(u)  is  continuous;  see  [l3]  ,  p.  25.  This  incidentally  shows 
that  G(0)  =  a.  If  V  3>  i  the  absolute  continuity  of  H(u)  follows  from 
the  Plancherel  theorem.  See  any  text  on  Fourier  transforms.  In  any 


‘m 


case,  define  the  functions  hm(u)  ■  — • L-  /  e“itu ^(it)dt,  m*l,2, 

Xni 

An  integration  by  parts^  gives  lor  u^O 
(4.1) 


h  (u)  =  — 

m  2TTiu 


'f(ira)e'lau-'f(-iT,)elraU  + 


»-  m 


e'itu  mtt)  dt. 
dt 


If  0<u1^u^u2,  (4.1),  (3-4),  and  (3,f)  show  that  the  continuous 
functions  hm(u)  converge  uniiormly  in  Ju^u^j  to  a  continuous  function 
h(u).  Moreover 

(4.2)  H(u2)-H(u1)  =  lira  f  (e  2-e~  vf/(lt)dt 

m — ^00  J_m  -/TTit 


the  first  equality  in  (4.2)  following  from  [13]  ,  p.  28  and  the  second 
from  the  fact  that  the  hm(v)  are  uniformly  bounded  for  u^<uS^U2« 

In  case  Ez^^ooand  r  <"/  f  k  -  1,  repeated  integration  by  parts  of 
(4.1)  ^nd  reference  to  remark  (b),  Theorem  (3.3),  sho»s  that  the  first 


I  am  indebted  to  J.  V.  Tukey  for  this  suggestion,  which  simplifies 
the  original  oroof. 
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r  derivatives  of  h(u)  are  continuous  if  u^-0.  The  usual  integral  ex¬ 
pression  for  h(u)  in  terras  of  Y^it)  shows  that  V  rfl  implies 
h^(u)  is  continuous  at  0. 

Corollary  to  the  continuity  of  H(u)  :  the  numbers  P(zn*r) — ^*0 

mllomh  ill  r,  r  >1,  ajS  n — *c£>-  We  have  pnr  =  * 

f(?)'°(?  -  pi)  ♦  C(s^)“°n(£j)  •  The  deslred  result  f0ll0W, 

because  G^(u)  y  G(u)  uniformly  for  u  0  and  because  G(u)  must  be 

uniformly  continuous  for  0  ^  u  <  CO  (right-continuity  at  0). 

We  next  consider  the  behavior  of  H(u)  near  u  —  0,  when  Theorem 

3.3  suggests  what  sort  of  result  may  he  expected.  If  the  function  M(s) 
of  Theorem  3.3  were  a  constant  M  it  would  follow  from  a  Tauberian 

y 

theorem  due  to  Karnmata  (see  [I4]  ,  pp.  139-192)  that  H(u)/~Mli _ 

r  cad 

as  u — *■  °r  p ?7»1)  ‘  In^e?ratl-nB  both  sides  of  this  re- 

lation  from  u  to  ux  would  give 


(4.3) 


Mdv 

v 


The  analogue  of  (4.3)  turns  out  to  be  true,  as  shown  by  Theorem  4«2, 
which  shows  that  in  a  certain  mean  sense,  H(u)  behaves  like  u^  U(u) 

rw~* *) 

as  u  - y  0t .  (Ve  defined  M(u)  =  V(-u)  for  u>  0.) 


Theorem  4.2.  Lim 

u  — y  0t 


ux 


u 


flixldY-  =  —4  r  M 


(v)dv 


The  proof,  which  follows  directly  along  the  lines  of  the  proof 
of  Kararaata»s  theorem,  is  sketched  briefly  in  section  9,  for  a  somewhat 


Revised  7-29-48 


-14- 


more  general  situation. 

A  corollary  of  Theorem  4.2  Is  that  if  V  <1,  h(u)  cannot  be 
bounded  as  u — wOfj  for  h(u)  <  K  implies 

— i —  /  Mill  dv  3>  0,  or  lim  u  !>0,  which 

u-t°+ju  77^1  ^  (7TT) Jx  v  u_ >0* 

implies  V>1.  An  example  to  be  given  in  section  5  shows  that  if 
V  :  1,  h(j)  is  at  least  in  certain  cases  bounded  but  discontinuous 
at  0. 

In  order  to  consider  the  behavior  of  H(u)  as  u — >-00  we  first 
prove  a  theorem  which  applies  to  any  distribution  whose  m.g.f.  is  an 
entire  function. 

Theorem  4-3*  Let  F(u)  be  any  c ,d.f .  whose  n.g.f ,  <^(s)  is.  entire. 
Let  p  be  the  order  of  £^(s).  Let  Q  be  defined  by 

q 

dP  (u)  <00. 

Then  .  t  _L_  =  1. 

The  proof  is  given  in  section  9. 


Q  *  l.u.b.  q : 


•00 


'-00 


u 


Combining  Theorems  3.4  and  4*3,  we  obtain  immediately 


Theorem  4.4.  Let  Q  l.u.b.  q: 


9 


f°A(u 

Jo 


)du  <0o*  Then 


Here  is  given  by  definition  3.4.  If  f(^)  is  not  a  polynomial, 
whether  entire  or  not,  the  proof  of  theorem  4.3  will  show  that  Q  =  1, 


Pefore  completing  the  present  proof,  the  writer  communicated  this 
result  to  P.  P.  Boas,  Jr.,  who  sent  back  a  proof  along  different 
lines. 
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_  k 


and  we  Interpret  theorem  4*4  in  that  sense.  The  trivial  case  f(s)  -  s 
is  excluded,  so  /°>i. 

Conjecture .  Let  £  ( s )  of  Theorem  4*3  he  of  finite  order  and 

r0o  .  ti  iQ 

type  C,  0  <C  <O0.  Let  Q  =  and  let  A  =  l.u.b.  A*  ij  e  dF(u)<f<>0. 
Then  (CP)Q  •  (AQ)  =  1. 

The  proof  for  the  case  rational  follows  the  same  lines  as  the 
proof  of  Theorem  4.3;  a  general  proof  has  not  been  found.  If  the  con¬ 
jecture  is  true  then  having  determined  and  Q,  v.hen  k(s)  is  a  poly¬ 
nomial,  and  having  estimated  C  by  the  procedure  indicated  following  the 
corollary  to  theorem  3.4,  ve  obtain 


(4.4) 


■  i 


TOO  a»uQ 

for  the  l.u.b.  of  the  numbers  A’  such  that  /  e  h(u)du  • 

Jo 

u 

The  corresponding  number  A  which  applies  to  g(u)  is  given  by 


(4.5) 


A*  =  A(l-a)  . 


5.  Some  special  cases .  In  this  section  ve  shall  discuss  some  special 
cases  in  which  the  m.c.f.  0(s)  and  the  c.d.f.  G(u)  may  be  determined 
explicitly.  For  these  cases  and  for  certain  others  there  is  a  close 
relationship  between  the  simple  discrete  branching  process  and  another 
type  of  model  to  be  discussed  in  section  8.  Finally  a  numerical  com¬ 
putation  of  the  distribution  G(u)  vill  be  given  for  a  particular  case 
where  f(s)  is  a  second  degree  polynomial. 

Suppose  f(s)  has  the  form 


f(s)  1 '  k +  a;(u<x-c?s) 
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2 

with  x>  1,  06  x  -  1,  where  f'(l)  =  x  end  f  (1)  ♦  f  (1)  2  Ez^  • 

x(l4-2CQ.  It  Is  easily  verified  (as  pointed  out  by  ?oincare/  in  }ll[ ), 
that  the  solution  of  the  equation  0(sx)  =  fp(s)]  is  given  by 

0(s)  -  1  +  y;ith  0(0)  =  0  (0)  *  1.  The  number  a  satisfying 

x-1-  03 

a  =  f(a)  is  given  by  a  =  C6^1 -X .  The  functions  s )  and  k(s)  of 

section  4  are  given  by  M^(s)  =  1 . ,  k(s)  =  _  s  .  The  number  7 

'  1-s  x-(x-l)s 

of  Theorem  3.3  is  1.  The  density  function  h(u)  (definition  4.1)  is 

simply  e~u,  as  seen  by  direct  calculation.  The  number  Q  of  Theorem 

4.3  is  1,  as  it  should  be,  since  f(s)  is  not  an  entire  function.  The 

c.d.f.  H(u)  is  l-e“u,  and  H(u)  /'-wu  near  u  =  0,  in  agreement  with 

Theorem  4.2.  Various  aspects  of  the  case  f(s)  =  As  +  _B  have  been  dis- 

CsTd 

cussed  by  numerous  authors. 

Somewhat  more  generally,  we  may  consider  generating  functions  of 
the  form 

(5.1)  k(s)  -  S[l-(»-l)s‘]'1/’1  X>1. 


The  function  k(s)  is  a  generating  function  if  and  only  if  m  is  a  non¬ 
negative  Integer.  In  this  case  we  have  0(s)  ■  v('/(s)  -  (1-ns)  and 


g(u)  =  h(u) 


Here 


end  we  note 


that  unless  m  = 1 


thedenslty  function  h(u)  is  unbounded  near  u*0.  A  physical  interpretation 
for  this  case  will  be  given  in  section  8. 

2 

As  a  numerical  illustration  we  consider  the  case  f(s)  *  0.4s  +  0.6s  . 

2  2 

We  have  x  r  F.z^  =  1.6  and  <J  *  E(z^-x)  =  O.24.  For  the  asymptotic 


distribution,  Fw  *  1, 


E(w-l)2  =  ^L 

X  -X 


0.25. 
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The  number  "/  -  log^  f  (0^4)  =  1*9495  so  that  ^(s)  which  is  identical 
with  0 ( s )  in  this  case,  isQ/ 


1.9495 


as  s  goes  to  Oo  with 


near 


Fe(s)^0.  This  Implies  that  the  c.d.f.  H(u)  and  likewise  G(u),  since 

r  .  ,  A~|  /  .  1.9495 

the  tv.o  are  equal  here,  behaves  like  1/p  (1+V)  M(u)  times  u 
u  =  0,  where  the  "behavior"  is  in  the  sense  of  Theorem  4.2.  Numerical 
determination  of  M(u)  would  not  be  difficult.  The  number  ^  of  Theorem 
4.4  is  given  by  log  2  -  I.4748.  This  means  that  tt*)  is  an  entire 

X 

function  of  order  I.4748  and  hence  that  the  density  function  h(u)  goes  to 


-u 


Q-e 


-u 


q+£ 


?ero  more  rapidly  than  e  and  less  rapidly  that  e  for  any 

6  >  0,  where  Q  =  =  3.1061,  and  "more  rapidly"  is  used  in  the  sense  of 

Theorem  4.4. 

The  function  L(s)  =  lim  1°&  (i*  ,).  was  computed  for  four  values 

n — hjo  s/°2n 

of  s  between  s  =  1  and  s  -  x  -  1.6;  in  each  case  the  value  was  0.744625 

so  that  it  appears  likely  that  in  this  case  L(s)  is  constant.  Hence 

C  =  Max  L(s)  =  0.744625  and  the  quantity  A  defined  by  (4.4)  is  0.26430. 

Thus  the  conjecture  following  theorem  4.4  indicates  that 

3.1061 


f 


g(u)  e 


(0.744625  -  £  )  u 


du  is  (divergent,  convergent)  according 


as  the  +  or  —  sign  holds. 

Through  the  kindness  of  Mr.  Cecil  Hastings  of  the  Douglas  Aircraft 
Company,  the  c.d.f.  G(u)  was  computed  for  this  case.  The  coefficients 
in  the  power  series  expansion  of  0(s)  were  obtained  from  the  functional 
equation  (3.1)  and  G(u)  was  then  obtained  by  inverting  0(it).  The  values 
of  G(u)  are  given  in  TabLe  I. 
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TABLE  i 


G(u), 


the  limiting  probability  that 


u  for  the  case 


f(s) 


0.4s 


♦  0.6s 


2 


u 

G(u) 

0.00 

.00000 

0.25 

.04753 

0.50 

.37275 

0.75 

.34550 

1.00 

.53117 

1.25 

.60932 

1.50 

.83042 

1.75 

.91357 

2.00 

.96781 

2.50 

.99751 

3.00 

.99993 
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6.  Number  of  generations  to  extinction.  It  v,as  pointed  out  in  section 

2  that  when  x  <  1  the  probability  is  1  that  zr  =  0  for  some  integer  n. 

We  assume  through-out  section  6  that  x  <  1. 

Definition  6.1.  Let  the  random  variable  N  be  the  smallest  integer 

n  ^ueh  that  z  ,  3  0.  Define  the  moment-generating  function  of  N  by 

ntl 

CD 

0(S)  --  X=  ensP(N-n). 
n=0 

00 

Clearly  P(Nm)  =  -  PnQ,  so  that  0(s)  *  gT-  e"8(Pn+1,0- P„0)- 

Definition  6.2.  Let  bn  =  1  -  Pnt]^0>  v.ith  bQ  =  1  -  pQ.  The 
numbers  bn  satisfy  the  recursive  relation 

(6,1)  bn+1  3  1  -  f[l-bn]. 

Define  the  function  0^(s)  by 

Oo 

0,  (s)  =  2IL  b  enS. 

1  n=0  n 

We  see  that 

(6»*)  ©(s)  =  1  t  (es-l)@1(s), 

so  that  it  suffices  to  determine  the  function  ©^(s). 

The  function  @^(s)  belongs  to  a  type  which  has  been  studied  by 

Fatcu  M  and  Lattes  jl6]  .  If  we  let  es  =  z  we  see  that  ©1(z)  is 

a  power  series  whose  coefficients  are  successive  iterates  of  the  function 
f*(h)  =  l-f(l-b);  i.e.,  b ^  =  f*(bn)  =  *here  f#(°)  =  °> 

#  t 

f  (0)  =  x  <^1.  It  vas  shown  by  Fatou  that  a  function  of  this  sort  is 


(Revised)  7/29/48 


-20- 


meromorphic  with  poles  at  s  r  -n  log  x,  n-l,2,...# 
©iCs)  in  the  form 

8,0 

L  s  .  ^  s  ,  3 


An  expansion  for 


3 


1-xe 


1— x  e' 


l-x3eS 


was  obtained  by  Lattes,  the  expansion  converging  everywhere  except  at 
the  poles.  The  quantities  an(*  yQ  are  defined  as  follows:  the 

function  ^(s)  *^3  s2  */4<s3  4  ...  is  determined  by  the 

X  2 

functional  equation  v^(sx)  =  f  [y^(s)J  with  the  condition  /*'  (l)  1, 

The  number  yQ  is  determined  by  ^(y  )  *  bQ  -  1  -  pQ.  Perhaps  the  easiest 

way  to  determine  y  is  to  use  the  fact  that  the  Inverse  function/^  x(s) 

o 

satisfies  the  functional  equation^ “1  [f*(s)]  =  xv^^s),  from  which 

we  can  determine  the  power  series  for  (bQ) . 

Since  the  use  of  Lattes*  sxpansion  requires  finding  the  expansions 

of/^(i)  An6y^~^(s),  we  now  give  another  method,  giving  a  different 
kind  of  expansion;  this  method  appears  particularly  adapted  to  the  ease 
here  illustrated,  where  f(s)  is  of  the  second  degree.  Then  (6.1)  becomes 

<6»3) _ brul_*  xbn  -  p2bn'  bo  "  1  ‘  V  ' 

PiflnlUgn  6.3.  The  functions  0  (s),  k=l,2,...,  are  given  by 

ML 

(6-4)  0k(,)-g(bn)ken3. 


If  we  raise  both  sides  of  (6.3)  to  the  k  power,  multiply  both  sides 
ns 

,  sum  on  n  irotn  u  to  ». 

k 


by  elia,  sum  on  n  from  0  to  Oo  ,  and  solve  for  0,  (s),  we  obtain 

•  1/  ' 


(6.5) 


ek(s) 


-  v"*SO) 


e"s  -  x* 
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(Justification  for  the  rearrangement  of  series  will  come  out  of  the 
subsequent  proof.)  If  we  put  k  =  1  in  (6.5)  we  obtain 

(6.6)  e.(s)  .  V  S  ~  p*92(s). 

1  -s 


(s) 

(s) 


Definition 


and  Fn(s), 

.  Vs 


6.4.  We  define  recursively  sequences  of  functions 
such  that  for  each  n,  ©^(s)  =  Sn(s)  ♦  ?n(s).  Let 


P1(s) 


P2e2(s) 


Suppose  now  that  Rn(s)  c**  th# 


form  An^®n+^(s)  t  ...  ♦  A  nn®2n^s^*  the  AnJ  being  factions  of  s,  p^, 

and  x,  but  not  explicitly  of  bQ;  while  Sn(s)  is  a  rational  function  of 
e-S,  P2,  and  x,  and  a  polynomial  of  degree  n  in  bQ.  Now  put  k  =  n  +  1 

in  (6.5)  and  substitute  the  expression  obtained  for  ©^^(s)  into 
Fn(s).  Collecting  terms  we  now  define  Rn4>1(s)  as  the  sum  of  terms 
involving  en+2(s) . 02n+2(s):  Hntl(»)  =  *ntl,i©nt2(»)*- 

An+I,n+I@2n+2(S)!  th€n  Sn.l(s)  =  0l(s)  '  Vl(s)  13  3  rational 
function  of  e”s,  p_,  and  x,  and  a  polynomial  of  degree  n  +  1  in  b  • 

2  O 

2 

Theorem  6.1.  Let  f(s)  =  Pq+P^s+p^s  »  vith  x  <  1*  Suppose  that 

x  +  p^b  <  1.  Then  the  functions  S_(s)  converge  to  0, (s)  in  a 
2  q  n  JL  — 

neighborhood  of  s  -  0. 

The  restriction  x  +  p2b0  <  1  may  fail  to  hold.  Hovever  thi3  is 

not  a  serious  restriction;  we  pick  a  value  of  n  so  that  x  ♦  p2b  <T 1. 

^  \  v.  (n-l)s  ns~*  v  n 

Then  ti^s)  =  bQ*. .  •♦bn_ie  ♦  e  ©x(s),  where  the 

00 

function  0.  (s)  =  2 _ b  n^S  is  the  same  type  of  function  as 

1  J=n  J 
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Q  (s)  $  theorem  6.1  Is  then  applicable  to0^(s). 

If  the  conditions  of  theorem  6.1  are  satisfied,  we  have 
(6.7)  ©1(s)  =  boe”S|pf1(s,x)  -  (s,x)  ♦  2xp2  bjV)^ (s>*) 

-  P23b03^e"bf  5x3)TT4(s,x)f ...J 
k  /  \  t 

where  TT  (s,x)  =  TT(  -g1-  — )  •  Since  E(N)  =  ©  (0)  =  0.(0)  ,.nd 
k  r-L\e-s.xT  I  1 

E(N2)  =  0"(O)  =  ?0{(O)  *  ©1(0),  we  have 

E(N)  =  bojTT1(o,x)  -  PkboTT2^0,x^  *  2xp  fbfofe.x) 


-  P23bo(le5x3)Tl^(0,x)t...']  , 


E(N2)  =  -  F(N)  ♦  2brTT[(0,x)  -  P^TT^O.x) 

t  2xp2b3rr(o,x)  -  (5x3+i)  p3b3rr(o,x ) 

2  °  3  c  o  4 

♦  p  3b 

2  o  4 

i  k 

where  TT,(0»x)  =  TT  (0,x)  >_  — —  . 

k  k  r=l  i_xr 

We  now  prove  that  if  x  +  p^b^  <T  exPansion  (6.7)  is  valid 

in  some  neighborhood  of  s  -  0.  We  shall  denote  the  particular  values 
of  x,  p^,  and  bQ  vith  which  we  are  dealing  by  x,  p^,  and  bQ.  Now  let 
x,  p  ,  and  b  .  be  three  complex  numbers,  arbitrary  except  for  the  fol¬ 
lowing  restrictions: 
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(6.8)  x  ♦  p  <  1,  b  <  1 

2  0 


and  define  the  numbers  b  In  terns  of  b  ,  x,  and  p  ,  by  means  of  (6.3), 

n  o  2 

vith  0  (s)  defined  by  (6.4). 

K 

%e  first  show  that  (6.7)  is  valid  if  (6.8)  holds,  and  then  show  that 

the  domain  of  validity  ilso  includes  the  original  numbers  x,  "p,  ,  and 

b  ,  provided  "x  +  pLF  <1. 
o  *  o 

If  (6.8)  is  satisfied,  we  have  J b^  A  x|  where  A  is  a  positive 
constant.  Nor  suppose  1  T  Then  the  series  defining 

©  (s),  k-1, ?,...,  are  uniformly  and  absolutely  convergent  in  the  domain 

J\ 

| e s |  <  T.  Moreover,  if  |x|  +  p^  -  A  <C  1,  ve  have  |bn  5ibQZ\n 
rhence,  if  k  is  an  integer  large  enough  so  that  TZ\  ^  , 

(6.9)  G.  (s)  <-2bk 

k  o 


IS 

e  <>  T.  Now  write  0^(s) 


A  .(p.  ,x,s)Q  .(s),  vhere  n  is  large  enough  so  that 


Tz^\n  <  -.  Let  A  (n  , x , s )  =  Max  L  .(p0,x,s).  Passing  to  the  next 
2  n  2 9  l^j^nl  nJ  2 

.  A  J  A  nfl  ^  (  /v  ntl  \ 

stage  ve  see  that  A  ,  <=  A  +  1  nil  2-\  S*  A  Ilf  — *  I  . 

n+1  n  -s  n  n  \  -s  n  / 

e  -x  \  e  -x  / 

Hence  the  numbers  A  are  bounded.  This  fact,  together  with  (6.9), 

n 


shows  that  lim  Fn(s)  3  0. 
n — >oq 
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Now  suppose  that  x  and  bQ  have  their  original  values  x  and 


while  p  is  small  enough  in  absolute  value  so  that  x  ♦ 


<  1.  In 


this  case  lira  Sn(s)  =  @^(s).  Ye  observe  that  S  (s)  is  a  polynomial 
n — y  oq 

of  degree  n  -  1  in  p^  and  that  S  ^(s)  is  obtained  from  Sn(s)  by  adding 
a  single  terra  of  degree  n  in  p„ .  Thus  ©^(s)  has  been  expressed  as  a 
power  series  in  p? .  Nov,  consider  0^(s)  as  a  function  of  p,  ,  with 


b  =  b  ,  x  —  x .  If  x  wE 
o  o  o 

is  analytic  in  p^  for 


<  1,  we  have  b^  =  0[(x)n]*  Thus  ®x(s) 

_= —  and  the  expansion  in  (6.7),  being  a 
o 


power  series  in  p„,  must  be  valid  when  x  t-  p,.  b  1. 

<-  cO 


7.  Estimation  of  parameters .  Until  nov.  we  have  assumed  that  the 
parameters  pr  are  known  numbers.  Ye  may  vish,  however,  to  estimate 
them,  having  observed  the  numbers  , . . . , z ^ .  In  order  to  get 

simple  maximum  likelihood  estimates  for  the  p^,  it  appears  necessary 

to  introduce  certain  auxiliary  random  variables. 

Definition  7.1.  Let  z  .  be  the  number  of  individuals  in  the 

mk 

generation  who  have  exactly  k  desoendents  in  the  (mfl)st  generation. 

Let  Z  =  ltz, t . . . +z  . 
n  1  n 

Theorem  7.1.  Maximum  likelihood  estimates  of  pr  and  x,  b^sed  on 
observed  values  of  z  .  for  ra  ^n,  are  respectively. 

- -  - - mlc  -  >  -  - L - - 


(Note  that  the  estimated  involves  only  z^, . . . • ) 
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If  z  is  fixed  the  Joint  conditional  probability  function  of 
m 

CO  7  -I  I  co 


z  ,z 
mo  m 


r-  w  z  *1  /  vs> 

. . ..,  is  (z  )!  TT  p  i  /TT  (z  )*  •  Thus  the  Joint  probability 

1  L  m  r=0  r  J/  r-0  mr 


function  of  the  z  for  ra=0,l, . . . ,n,  and  r=0,l,2,...,  is  given  by  the 

product  of  two  factors,  one  of  which  is  independent  of  the  p  ,  the 

to  /  n  \  r 

logarithm  of  the  other  being  z (5  z  ]  log  p  .  The  value  of  this 

r=0  V^O  mTJ  r 

/\ 

expression  is  clearly  maximized  by  taking  pr  =  pr  as  given  above. 


Since  V  z  -  z  and  S  rz  =  z  ,,  the  quantity  >  rp_  gives 
^ —  mr  m  ^ —  rar  mtl  —  r 


/r  ns  above. 


Although  the  estimates  p  are  the  same  as  we  would  obtain  if  we 

r 

were  dealing  with  Z  trials  from  a  multinomial  distribution  with  pro- 

n  n_ 
babilities'  p  ,  the  Joint  distribution  of  the  quantities  ^T_  zmr> 


m=0 


r-0,1,...,  is^not  multinomial.  nFor  example,  if  Zn  >  1  the  probability 


r  n 

of  the  event m  z  -  Z„ ,  ^  _  z  =0  for  r 

l^m-0  mo  n  m=0  rar 


r^Q J  is  0 


We  shall  next  show  that  the  estimate  x  is,  in  a  certain  sense, 
consistent . 

Theorem  7.2.  If  x  1 ,  the  random  variables  Z^^  J converKe  in 

#  *  i  * 

probability  to  the  random  variable  xV  where  V  =  i  if  w  *  0  and  V  »  1 


if 


0. 


If  w^.  0  then  for  all  n,z  =JL  0  and  1  Z  — >-Ooas  n - WCO» 

7  n  /  n 

Hence  in  this  case  ^7^ ^-ljy/zn  converges  to  x  if  Z^^y/zn  does.  On 

the  other  hand,  P(w  =  0)  -  a  -  P(zn=0)  for  some  n,  so  that  if 
w  =  0,  ZfM.1  /Zn  =  1  with  probability  1  for  n  large  enough.  Thus  we 
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need  only  show  that  converges  to  x  if  x  >1  and  *  ^0. 

We  need  the  foll^.ing: 

Lemma  7.1.  if  x  >  1,  the  random  variables  xn  converge  in  pro¬ 


bability  to  wx  . 

x  - 1 


/ 


Since 


(7.1) 


wx  - 


x-1 


n 


n 

n 


w 

1TTT 


— r  - 

x - 1 /  r=0 


(w-wp) 


,n-r 


it  will  be  sufficient  to  show  that  lim  /  x 

n  —  >oq  \x-l; 


1  E(w^)  =  0  and 

2nt2 


lira 


n — ►  © o  \  r=0  xn_r 


(v  “V  ) 

'  =  0.  The  truth  of  the  first  statement  is  obvious, 


since  Ew*  is  finite.  It  follows  from  (2.5)  that  E(w  w  )  -  Ew  if  s  >  r, 

r  s  r 


r-~ 

E(ww  )  =  lim  E(w  w  )  =  whence  E (w-w  ) 


I2L 


n — foo 


(x2-x)xr 


and 


2 

Ef(w-W  )(w-w  )]  -  __£L? -  if  s  >  r .  Then  E  ^W~WrA  - 

L  r  s  -J  tJ*  l  -n~r  I 


(x  -x)x 


_l_  sfL 

2n 

X  x  -X 


"  n  n  s-1 

z:  *r  *  22Z 

r-0  s - i  r=0 


,  and  this  quantity  clearly 


approaches  0  as  n — V  CC»,  proving  Lemma  7.1. 

Define  the  random  variables  w  and  V  as 

n 

w  =  v  when  v.  =£  0 

w*  =  L  when  w  =  0 

V  =  v  hen  z  =k  0 

n  xn  n-T 


n  7.l 


x  v  hen  z„  —  0 


n 
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It  is  clear  that  the  V  converge  in  probability  to  w  -JL.,  and 

n  x-1 

#  # 

we  note  that  the  c.d.f.  of  w  is  continuous  at  w  =0.  Hence, 


Lim  P 
n — >Co 


>  e  >  o 


lim  P 
n — VOO 


(v  ,-v  ±v„£  >  o  ]  = 

\  n+1  n  n  ^  ) 

p(  - JL  *£  ^  0]  =  o.  It  follows,  under  the  conditional  hypothesis 

\  x-1  >  / 

w^O,  that  the  variates  ^n-fl  converge  in  probability  to  x,  since 


n 


"ntl  =  x  ^ntl  when  zn^  ^  0* 


n 


n 


8.  Continuous  Models.  As  mentioned  in  section  1  there  are  situations 
where  it  is  more  Important  to  consider  the  number  of  Individuals  existing 
at  a  given  time  than  the  number  in  a  given  generation.  Let  a  set  of 
probabilities  pr  be  given.  The  question  arises  whether  we  can  interpret 
these  as  probabilities  that  an  individual  will  have  a  given  number  of 
descendents  at  the  end  of  some  fixed  period  of  time.  V.’e  might  then 
suppose  that  each  individual  in  existence  at  that  time  has  the  same 
probabilities  of  having  a  given  number  of  descendents  at  the  end  of  the 
next  (equal)  length  of  time,  these  probabilities  being  independent  of 
the  age  of  the  individual.  A  model  of  this  sort  might  be  considered  in 
certain  fission  processes,  if  the  probability  of  fission  is  independent 
of  age.  It  should  be  noted  that  the  "descendants"  of  an  individual  may 
Include  the  individual.  For  example,  if  a  bacterium  splits  in  two  we 
may  either  regard  it  as  having  produced  two  descendents  and  dying,  or 
as  having  produced  one  descendent  and  itself  surviving. 
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If  an  interpretation  of  this  sort  is  to  be  satisfactory,  inter¬ 
polation  in  time  must  be  possible.  In  other  words  there  should  exist 

i 

a  family  of  functions  f^(s)  defined  for  all  positive  n  such  that 
fn  =  *n  +n  such  that  for  each  positive  n,  fn(s)  is  a 

C0_ 

probability  generating  function,  fn(s)  =  ]>_  P  (n)sr;  and  such  that 

r=0 

for  n  =  0,1,2,...  the  functions  fn(s)  coincide  v.lth  the  iterates 

, . . .,  V’e  nay  then  interpret  fn(s)  as  the  generating 

function  at  time  n.  It  is  readily  seen  that  in  general  such  a  family 

of  functions  will  not  exist.  For  example,  if  such  a  family  exists  we 

must  have  f(s)  -  n^  iterate  of  f  .  (s)  for  arbitrarily  large  integral  n, 

1/n 


2. 


so  that  f(s)  cannot  be  a  polynomial  of  degree 

The  functional  equation  0(sx)  =  f|0(s)]  shows  that  f  (s)  =  0 ^x0-i (s)j  , 


Pll  J  u 

whence  fn(s)  =  0  xn0  (s)  for  integral  n.  The  expression  0|xn0-i(s)| 


then  can  be  taken  as  the  definition  of  f  (s)  for  all  positive  n.  See 

n 


Hadamard,  |V]  .  The  problem  of  determining  whether  the  functions  so  de¬ 
fined  are  a  family  of  generating  functions  will  be  discussed  in  a  sub¬ 


sequent  paper.  We  remark,  however,  that  if  f(s)  has  the  form  s 

x-(x-l)s 


considered  in  section  5  then  the  iterates  f  (s)  have  the  form 

n 

s _ ;  they  are  clearly  generating  functions  lor  all  positive  n, 

xn-(xn-l)s 

satisfying  the  required  relation  f  (f  )  =  f  .  Now  suppose  g(s) 

is  some  function  such  that  the  function  f(s)  =  Z 
generating  function  for  all  x^>  1,  with  g(l)  =  1.  As  pointed  out  by 
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Ulam  and  Hawkins,  the  iterates  of  functions  f(s)  of  this  form  are  con¬ 


venient  to  vork  vith,  the  nL  iterate  being  simply  g 


_ ti 

xn-(xn- 


d«(i 


In  addition,  the  requirement  that  f(s)  be  a  generating  function  for 

all  x  ■>  1  shows  that  the  functions  f  (s)  are  generating  functions  for 

n 

all  n>0.  The  simplest  function  g(s)  which  satisfies  our  requirements 
is  g(s)  -  sm,  where  ra  is  any  positive  integer.  In  this  case  f(s)  has 
the  form  considered  in  (c.l)  and  fn(s)  =  s  (xn-l)  s^J  ^  .  As  n - >0 


ie  have  fn(s)  =  (1  -  jj  log  x)s  + 


n-ffi  x  sm *l  t0(n2). 


We  may  inter¬ 


pret  this  as  follows.  A  particle  in  existence  at  a  given  time  may,  in 
a  short  time  interval  ^  t,  either  split  into  nul  particles,  with  pro¬ 
bability  -At_lo£_x;  or  it  may  remain  unaltered,  with  probability 

ra 

1  _  _^_t_JjOiL_Z  #  if  it  splits,  each  particle  produced  has  the  same 

chances  for  splitting  as  its  parent,  etc.  Thus,  from  the  results  of 
section  5,  it  follows  that  if  we  begin  with  a  single  particle  at  time 
t  -  0,  the  asymptotic  probability  density  function  for  where 

z  is  the  number  of  particles  at  time  t,  is  given  by 

w 


(m^V^)/r(*) 


It  is,  of  course,  customary  to  begin  with  the  elementary  probabi¬ 
lities  for  a  certain  number  of  births  in  a  short  time  At  and  determine 
the  functions  f n ( s )  from  these  Dy  means  of  differential  equations.  See, 
for  example,  Arlev,  JL7]  .  The  results  of  the  present  paper  can  be  ap¬ 
plied  in  some  cases  to  the  continuous  problem  even  when  an  explicit 

determination  of  the  f  (s)  is  difficult.  A  discussion  will  be  given  in 

n 

a  later  paper. 
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9.  Cone  proofs.  We  give  in  this  section  proofs  for  (A)  theorem  3.3, 

(B)  theorem  3.4,  (C)  theorem  4**,  and  (D)  theorem  4.3;  in  certain  cases 
re  shall  indicate  slightly  more  general  results. 

(*\)  We  make  use  of  a  result  of  Koenigs,  in  the  form  applicable  here. 

Koenigs’  theorem:  If  |s|  ^  A  1  and  then 

^(s)  =  q^  B(  s  )  |7*fQ(q^)J  where  B(s)  analytic  for  |s|  <^Aand 

satisfies  the  f unctlonal  e lua tion  E Ik ( s )J  -  q  t(s). 

I—  1 

Here,  Q(q^) 

We  remark  that  P(s)  ^  0.  The  proof  of  Koenigs’  theorem  follows  readily 

if  re  write  k  (s)  =  q^_ik(s)  TTu  +  .f^L  ,  where  £(s)  =  -  Q.  • 

n  1  !  '  n  si 


means  bounded  by  Aq^,  where  A  is  Independent  of  s. 


j=l|  "1 


Nor 


iv  let  t^  be  a  positive  number  such  that  j  Y(s)|  <  1  when 
0  <C  s|  ^  t  and  Be(s)  ^0.  (For  the  rest  ol  this  proof  re  as 


sume 


Fe(s)  0. )  Such  a  number  exists;  on  the  imaginary  axis  re  have 
Y(lO  “It  it  -  iE^r’)  Jt^  4  o(t2)  v.here  E^r' )  1,  w’  having  the  dis¬ 
tribution  branching  from  k(s),  shoring  that  |  <C  1  if  t^0  and 

sufficiently  small;  while  if  Fe(s)  <Co  re  refer  to  the  expression 


rOo 


y(s)  -  J  e  udH(u).  bPt  A  =  M fi x  |  ( s ) |  1  or  t ^j: 


ij  i  y  i 

>  t.  let  N(s)  be  the  smallest  integer  such  that  s  fx  ^  t 


V 

N(s) 


Then  vg(s)  -  It 


N(s) 


bH[,  tO(qi(=)^]-  Now  B[r(=»j]  r  qiBYs)J- 

Let  M(s)  =  s|^  B  ^Y(s)|  •  Then  M(sx)  =  M(s).  Also  log^ 


s  t, 


1  t  log. 


5  t^  ,  and  theorem  3.3  follows.  Clearly  M(s)  s  i 


<  N(s)< 
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continuous  for  t  x  ^^|s|  ^  t  and  hence,  by  functional  continuation, 
wherever  Pe(s)5^0,  s=^4o. 

Concerning  the  remarks  following  Theorem  3.3  we  have  the  following: 

(a)  If  Ez^  Co ,  r-fold  differentiation  of  ^(sx  n)  =  *n  ^ 

gives,  for  jsj.^  >  0, 


CM) 


where  Qr^  is  a  polynomial  in  •  Now 


kn(s) 


be  true  of 


=  0  (q^)  when  |s|  <  A;  because  of  analytic! ty,  the  same  must 

(J) 


knv'(S) 


defined  above.  Since 


Put  n  -  N(s)  in  (9.1),  N(s)  being  the  integer 
kN^  sOh[)  z  O^y^js  "*)j,  remark  (a) 


follows . 

(b)  B(s)  is  clearly  ->0  when  s>0;  hence  M(s)S>:0  when  s<0. 

Since  B(0)  =  0,  B(s)^-0  for  sufficiently  small  s^40;  since  ^(s) - VO 

as  s| — Voq>  M(s)^0  for  s  sufficiently  large;  since  M(sx)  *  M(s), 
remark  (b)  follows. 

(c)  If  V  -  CO  ,  i.e.,  q  =  0,  then  kn(s)  goes  to  zero  with 

great  rapidity  as  n — *-00,  if  s|  <C  1.  The  general  line  of  argument 
is  clear. 

(P)  Let  k(s)  be  a  polynomial  of  degree  d  >  1  with  real  coefficients 
k(s)  r  q^4- . . .  ,  with  a  non-negative  double  point ,  k (OQ  =  OL> 0,  and 

such  that  k(s)  C>  s  when  s  >  Ot.  Let  rco  be  any  solution  of  the 
functional  equation  'f'(ms)  =  k  [V(s)j  *hlch  is  contiruous  for  s  0 
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and  satisfies  'f'C s )  >  QC  tor  s  0 ;  here  m  Ajs  any  number  1 .  Then 
theorem  3.4  holds .  with  x  replaced  by  m. 

It  Is  not  difficult  to  shov. 


:*  that  if  0C  <C  s,  s  s~ ,  lim  k.(s) 

1  *  J — >-oo  J 

OO  uniformly  in  s.  Hence  ^(s) — ^  00  as  s  - ►•OO.  Write  F(s)  = 

d  -J 


log 


(1  +  1  y~  q  s’*5).  Then  d  niogY^(smn)  =  d  niog  kn[V(s)]  = 
*d  J  =  1 


taken 


(l-d"n)log  qd  /(d-1)  +  log  y(s)  +  ZZ  d' ,  s  being 

large  enough  so  that  F  (vi&H)  is  continuous.  Thus,  since  the 
functions  P(k^  ^  )  are  bounded,  the  functions  d  nlog  Kf^(smn)  con¬ 


verge  uniformly,  for  s  sufficiently  large,  to  a  continuous  function 


L  (s)  satisfying  L  (ms)  =  dL  (s).  Let  L(s)  =  t  /^L*(s),  where  ^ -  log  d. 

m 


m 


Theorem  3.4  now  follows  by  an  argument  similar  to  that  used  to  conclude 

OO 


theorem  3.3.  (Note  that  ^ _  d  ^  ^|^(s)j^  =  0(d~n)). 


(C)  In  order  to  avoid  negative  signs  we  work  with  the  Laplace  transform 
instead  of  the  m.g.f. 

Let  H(u)  be  nondecreasing  on  (0,00)  v 1 th  H(0)  =  0;  let  T(s)  -- 
■00 


ruu 

e~sudH(u)  be  finite  for  s  >  0.  Suppose  ^^(s)  =  M (s^  4  o(  -^)as  s — 


where  0  <C  V  <Tca,M(s)  i_s  continuous  and  satisfies  M(sx)  -  M(s)  for 


s>  o,  y  being  some  number  1.  Then  lim 


u- 


dv  =• 


p(V 


41) 


Mill  dv. 

V 
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Followinp  the  Lines  of  the  proof  of  Karamata's  theorem,  me  se® 

7-xy  y_ j_ 

s  ^H[s)ds  =  D  t  o(l)  as  s — ►‘CO  where 

y 

rx  rx y  y-i  r03  _su 

D  =  Mill  ds;  i.e.,  /  s  ds  /  e  dH(u)  =  D  ♦  o(l),  or 

J 1  s  >y  Vo 


-  su 

e  dH(u)  =  D  ♦  o (l) ,  or 


-su  -nsu 
e  e  dH 


replacing  s  by  (n+L)s ,  J~  ^  ^ds  j~  e”SUe  nSUdH(u)  r  D/(n*l)V+  o  (l)  * 
■ -b-S  —  f  e  Se  nSsV  ^ds  ♦  o(l).  It  follows  as  in  [1 4]  ,  pp.  189-192,  that 

I  cy)  j 0 

if  F(u)  is  any  function  of  bounded  variation  in  (0,1)  we  have 


(9.?) 


rxy  , 

lim  /  s' 

y - ooly 


i-r 

J  n 


-su  ,  - su.  .  . 
e  F(e  )dH(u)  = 


-SwV-1, 


e'°F(s  )sf  Ads. 


Let  F(e~s)  =  es  if  0 s  ^  1  and  0  otherwise.  Then  the  theorem  follows 
f rom  (9 . 2) . 

(D)  Theorem  4.3  is  true  if  F(u)  is  any  bounded  monotone  increasing 
function.  For  simplicity  ve  assume  that  F(l)  =  0;  it  is  readily  seen 
that  this  causes  no  loss  in  generality.  The  proof  is  given  for  the 
case  1  OO  >*  it  will  be  clear  that  P  -  1  implies  Q  =  OO,  while 

if  z9  =00  (or  if  £  (s)  is  not  entire)  Q  *  1. 

ruppose  m  and  n  are  positive  integers  such  that  m^/n  ^/^/(/^-l) . 

Then 


(°.3)  /  cxp(un)  dF(u)  = 


f,  /  u^dFtu)  <  n  C 


'rf  limit 


r=0  * 


f^O  (rn)!  (r-fl)a 
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where  c  =  -a - L\  interchange  of  integration  and  summation  are  Justi 

*  k! 


fied  by  the  positiveness  of  ail  terns  involved.  Suppose 
0  <  £  <  JL  -  (l-lj.);  for  k  sufficiently  large  the  inequality 

c  k  ^  ^  is  satisfied;  see  Id],  p.  253.  Hence,  using  Stirling’s 

formula,  we  see  that  the  last  series  in  (9.3)  is  dominated  by  a  series 
whose  rth  term,  for  r  sufficiently  large,  is  controlled  by  the  factor 

rrm^  ^  m  .  Since  l-^  +  £--§-  Is  negative,  the  series,  and  hence 
the  integral,  converges.  V'e  have  thu3  proved  ^  1 . 

Conversely,  suppose Let  J^(s  )  =  £  £k(s),  vhere 


cir±rns^trn  >  ^=0 , 1 , . . . , rn- 1 .  At  least  one  of  the  functions  <^k(s) 
r=0  trn 


k*0 


mu 


st  be  of  order  P  .  Ve  suppose  that  (s)  is;  if  not  the  argument 
would  need  only  slight  modifications.  T.’e  have 


°C30  m  (rm)  Ic 

(9.4)  /  exp(u~n')dF(u)  >  n  S  7=— — -1®- 

Jl  r=0  [(rtl)n]  l 


Suppose  0  <  £  <  i  -  •  From  [l8]  ,  p.253,  the  inequality 

rm  (  ^  +  C\ 

crm  ^  (rm)  ^  must  hold  for  infinitely  many  values  of  r.  As  in 
the  first  half  of  the  proof  this  shows  that  the  series  and  the  integral 
in  (9.4)  diverge.  Thus  1  >  1  and  the  proof  is  complete. 

If  /°  is  rational,  the  conjecture  following  theorem  4.3  can  be 
proved  in  a  similar  manner  making  use  of  a  relation  between  the  class 
of  an  entire  function  and  the  coefficients  of  its  series  expansion; 

see  I4J  »  P*  95 . 
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